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Abstract 


The NASA Nebraska Space Grant (NSGC) & EPSCoR programs have continued their effort to 
support outstanding research endeavors by funding the Numerical Simulation of the Combustion 
of Fuel Droplets study at the University of Nebraska at Lincoln (UNL). This team of researchers 
has developed a transient numerical model to study the combustion of suspended and moving 
droplets. The engines that propel missiles, jets, and many other devices are dependent upon 
combustion. Therefore, data concerning the combustion of fuel droplets is of immediate 
relevance to aviation and aeronautical personnel, especially those involved in flight operations. 
The experiments being conducted by Dr. Gogos’ and Dr. Nayagam’s research teams, allow 
investigators to gather data for comparison with theoretical predictions of burning rates, flame 
structures, and extinction conditions. “The consequent improved fundamental understanding of 
droplet combustion may contribute to the clean and safe utilization of fossil fuels” (Williams, 
Dryer, Haggard & Nayagam, 1997, U 2). The present state of knowledge on convective 
extinction of fuel droplets derives from experiments conducted under normal gravity conditions. 
However, any data obtained with suspended droplets under normal gravity are grossly affected 
by gravity. The need to obtain experimental data under microgravity conditions is therefore well 
justified and addresses one of the goals of NASA’s Human Exploration and Development of 
Space (HEDS) microgravity combustion experiment. 
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Rationale 


“The engines that propel missiles, jets, and many other devices are dependent upon 
combustion. Liquid fuel is sprayed into an engine chamber where it evaporates and bums, 
generating the thrust that propels the object forward” (Mashavek, 2001, ^ 1). The amount of 
thrust created depends on many factors, including pressure, temperature, the fuel droplet 
evaporation rate, and turbulence. Therefore, data concerning the combustion of fuel droplets is 
of immediate relevance to aviation and aeronautical personnel, especially those involved in flight 
operations. 

“The combustion of fuel droplets is an important part of many operations, such as the 
heating of furnaces for materials processing or home heating, power production by gas turbines, 
and combustion of gasoline in a car's engine” (Williams, Dryer, Haggard & Nayagam, 1997, If 
2). The Earth’s gravity prevents many theoretical predictions involving fuel droplet combustion. 
Additionally, drop towers and aircraft are unsuitable for this type of experimentation due to time 
constraints and unacceptable levels of microgravity. The experiments being conducted by Dr. 
Gogos' and Dr. Nayagam’s research teams, allow investigators to gather data for comparison 
with theoretical predictions of burning rates, flame structures, and extinction conditions. “The 
consequent improved fundamental understanding of droplet combustion may contribute to the 
clean and safe utilization of fossil fuels” (Williams, Dryer, Haggard & Nayagam, 1997, U 2). 

The present state of knowledge on convective extinction of fuel droplets derives from 
experiments conducted under normal gravity conditions. “Due to the increase in the extinction 
velocity with droplet diameter, under extinction conditions natural convection becomes 
negligible at large ‘droplet’ (porous sphere) diameters and important at smaller droplet 
diameters” (Bowen, Woods, Narayanan, Smith, & Gogos, 2000, 4.3.3 p. 1). As a result, any data 
obtained with suspended droplets under normal gravity are grossly affected by gravity. The need 
to obtain experimental data under microgravity conditions is therefore well justified and 
addresses one of the goals of NASA’s Human Exploration and Development of Space (HEDS) 
microgravity combustion experiment. 

NSGC & EPSCoR Background and Research Involvement 

The Nebraska Space Grant Consortium at the University of Nebraska at Omaha develops 
research infrastructure and enhances the quality of aerospace research and education throughout 
the state. This grant provides national leadership in applied aspects of aeronautics and allows 
Nebraska colleges and universities to implement a balanced program of research, education, and 
public service programs related to aeronautics, space science, and technology. The grant 
administers funds to recruit and train professionals for careers in the aerospace industry. 
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EPSCoR (Experimental Program to Stimulate Competitive Research) assists states with 
low levels of federal research and development support, thus responding to a Congressional 
concern about increasing the geographic base of federal research support. Nebraska EPSCoR is 
a statewide effort, which provides leadership for development of research and development in 
science and engineering throughout the state. Specific to the University of Nebraska at Omaha is 
the Aeronautics Education, Research, and Industry Alliance (AERIAL), a comprehensive, multi- 
faceted, 5 year NASA EPSCoR 2000 initiative. This contributes substantially to the strategic 
research and technology priorities of NASA while intensifying Nebraska’s rapidly growing 
aeronautics research and development endeavors. 

The partnership between the NASA Nebraska Space Grant (NSGC) & EPSCoR programs 
allows for the selection of outstanding research projects that positively impact aeronautical 
technology advancement. These programs have continued their effort to support such research 
endeavors by funding the Numerical Simulation of the Combustion of Fuel Droplets study at the 
University of Nebraska at Lincoln (UNL). 

The CEFD Concept 

Dr. Vehda Nayagam guides the Micro gravity Combustion Science Program in San 
Diego, CA. This program is conducting a project flight definition experiment to obtain data 
under microgravity conditions. The UNO CEFD collaborative research team (CRT) is 
developing a new comprehensive numerical model for the convective extinction of fuel droplets 
to validate this model. The data collected from each institution contributes to one of the long- 
term goals of the HEDS microgravity combustion program. Specifically, that which promotes 
“understanding that will permit lessons learned in microgravity combustion experiments and 
modeling to be used in optimizing combustion devices here on Earth.” 

A team of researchers from the University of Nebraska - Lincoln, led by Dr. George 
Gogos, is conducting a comprehensive computational study of fuel droplet combustion at 
atmospheric pressure and zero-gravity ambient conditions under forced convection. Through a 
collaborative effort with NASA Glenn Research Center, Dr. Gogos and his colleagues are 
developing a science education component that demonstrates how the combustion process 
changes due to microgravity. 

Simplified as well as detailed chemical kinetics are employed in the research. The 
studies provide insights that can be applied to improve liquid fuel combustion with greater 
efficiency and safety, and reduce environmentally-adverse effects. In view of the detailed 
chemical kinetics, substantial complexities and uncertainties are involved in modeling 
combustion of a moving droplet through the currently finding experimental research on 
combustion of a moving droplet through the Microgravity Combustion Science Program. 

The research focuses on the validated modeling of two key topics: a) Transient 
combustion of a moving droplet with simplified chemical kinetics; and b) Transient combustion 
of a moving droplet with detailed chemical kinetics. The first topic is currently being addressed, 
whereas the second one is a longer-term research project. This work is a direct extension of 
research funded by the NASA Nebraska Space Grant and EPSCoR Programs. “Dr. Gogos’ 
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studies on droplet combustion at atmospheric pressure include combustion of moving droplets 
with infinitely fast kinetics as well as with one-step global kinetics” (Bowen, Holmes, et al. 

1999, p. 19). 

Research success depends on the team’s considerable experience combined with recently 
published studies on comprehensive modeling of hydrocarbon oxidation, which employ detailed 
chemical kinetics. Dr. Gogos’ assistant, Dr. Daniel Pope, is also supported under the NASA 
Nebraska Space Grant and EPSCoR Programs. He has been working for over two years 
simulating combustion of a moving droplet with one-step kinetics and contributes immensely to 
the timely completion of the proposed work. 

Research Progress 

The current CEFD CRT research has focused on the development of a validated 
numerical model for droplet combustion in a forced convection environment. Funded by the 
previous five-year NASA EPSCoR grant, a quasi-steady numerical model, which utilized one- 
step overall chemical kinetics and a single diffusion coefficient to describe the mass diffusion, 
was developed to study the convective extinction of fuel droplets under zero-gravity conditions 
(Gogos, Pope, & Lu, 2001, p. 2). As a result of suggestions made in the review of the 2001 Pope 
and Gogos article and as a prelude to incorporating multi-step chemical kinetics schemes, the 
quasi-steady code is currently being modified to allow for the different mass diffusion 
coefficients associated with each pair of chemical species. This modification to the quasi-steady 
code is part of the systematic addition of modeling complexities that was presented in the 
original research proposal. The end goal of the research is to develop an experimentally 
validated droplet combustion model that can be used for accurate predictions of single droplet 
behavior in practical combustion systems. 

The conditions present in convective droplet combustion experiments are different from 
those present in practical combustion systems. Droplet combustion experiments under forced 
convection are conducted by suspending the fuel droplet from a silica fiber in an ambient 
oxidizer at a fixed temperature (Too) and pressure (poo), as shown in Figure 1 . The oxidizer is 
"blown" over the droplet at some fixed velocity (Uoo). If the ambient temperature is high 
enough, or if an external ignition source is present, the droplet will ignite. The initial flame 
configuration (wake, transition, or envelope) depends on the "blowing" velocity. In practical 
combustion systems, droplets are injected into a combustion chamber. This situation is shown in 
Figure 2, where the droplet is injected, at some initial velocity, into a stagnant environment at a 
specified temperature and pressure. The moving droplet experiences a drag force that opposes 
its motion and the droplet velocity decreases. The initial velocity determines the initial flame 
configuration. If the initial flame configuration is a wake flame, the decrease in droplet velocity 
can result in a change from a wake to a transition flame, and finally to an envelope flame. The 
numerical model must be able to deal with both the suspended droplet case (for model 
validation) and the moving droplet case (for practical predictions). A transient code is required 
to model the change in droplet diameter caused by evaporation at the droplet surface, the change 
in flame position and configuration, the internal heating of the droplet, and the decrease in 
droplet velocity for the moving droplet. 
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A transient code has been developed to model droplet combustion in a forced convection, 
zero-gravity environment. One-step overall chemical kinetics and a single diffusion coefficient 
to describe the mass diffusion were used in the model. The model has been validated using the 
numerical results of the 2001 Gogos and Zhang research for the evaporation (no combustion) of 
n-heptane droplets in nitrogen at atmospheric pressure. Excellent quantitative agreement was 
observed for various ambient temperatures and initial droplet velocities. 

The transient code has been used to numerically investigate the combustion of n-heptane 
droplets in air at atmospheric pressure. The initial droplet diameter (DO = 0.5mm), initial droplet 
temperature (TO = 297K), and the ambient temperature (Too = 1000K) were fixed and the initial 
droplet velocity (Uoo(0)) or "blowing" velocity (Uoo) was varied. Results have been obtained for 
moving and suspended droplets with initial Reynolds numbers of 10, 25, 50, and 100. At a given 
initial Reynolds number, the fixed "blowing" velocity (suspended droplet) and the initial droplet 
velocity (moving droplet) are equal. The numerical results indicate that the initial Reynolds 
number determines the flame configuration that forms during droplet ignition for both moving 
and suspended droplets. An envelope flame is formed during droplet ignition for an initial 
Reynolds number of 10 and a wake flame is observed at the higher initial Reynolds numbers. 

This is in qualitative agreement with the quasi-steady code, which predicts an envelope flame for 
Reynolds numbers less than 12 under these same conditions. Once the initial flame 
configuration had formed (either envelope or wake), the suspended droplet cases exhibited the 
same flame configuration throughout the droplet lifetime. In the moving droplet cases, the wake 
flame that formed at the higher initial Reynolds numbers, gradually approached, and then 
eventually surrounded the droplet in an envelope flame configuration as the droplet velocity 
decreased. The predictions indicate a marked difference between the behavior of suspended and 
moving droplets. 

The development of a transient droplet combustion code represents a significant step in 
our current research which is funded by the new five-year NASA EPSCoR grant. The 
modification of the quasi-steady code to include multiple diffusion coefficients is nearing 
completion. Once this modification is tested, it will be incorporated in the transient model. The 
next step will then involve the incorporation of multi-step chemical kinetics in the quasi-steady 
and transient models. 

Additionally, a computer code has been developed to study the combustion of liquid fuel 
droplets in a combustion chamber. Figures 3-5 show numerical results for the temperature 
distribution surrounding a n-heptane fuel droplet that is moving through a combustion chamber 
at three different velocities. The droplet diameter is 500 microns and the temperature and 
pressure of the combustion chamber are 1200 Kelvin and 1 atmosphere respectively. The color 
bar to the right of each figure indicates the temperature range (500 to 2500 Kelvin). Figure 3 
corresponds to the lowest velocity shown. At lower velocities, the flame (high temperature 
region) surrounds the droplet in what is called an envelope flame. As the velocity is increased 
(Figure 4), the flame extinguishes near the front of the droplet. Further increases in the velocity 
(Figure 5) cause the flame to move toward the rear of the droplet forming a wake flame. 
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2001-2002 Progress 


The CEFD CRT has made substantial progress in 2001. Two major aspects have been 
added to the validated axisymmetric quasi-steady code for the combustion of a moving droplet. 
First, the assumption of a single mass diffusion coefficient has been relaxed. The current code 
allows for different multicomponent diffusion coefficients. Second, the quasisteady assumption 
has been relaxed. The current transient code can predict, for example, how a moving droplet 
with a wake flame configuration, transitions to a droplet surrounded by an envelope flame, as the 
droplet decelerates due to the drag it experiences. These added aspects are presented in more 
detail below. 

I. Numerical Simulation of Droplet Combustion in a Forced Convection Environment; 

Multicomponent Diffusion Coefficients. 

Over the past few years, a quasi-steady numerical model has been developed to 
investigate fuel droplet combustion under forced convection in a low-pressure zero gravity 
environment. Two key assumptions used in the quasi-steady model were; (1) a one-step overall 
reaction was used to describe combustion, and (2) all binary diffusion coefficients were assumed 
to be equal. In the original code, the values for the activation energy, and oxygen and fuel 
concentration exponents in the finite-rate kinetics were adopted from Westbrook and Dryer. The 
pre-exponential factor was determined by comparison of numerical results with experimental 
data for extinction velocity. 

For n-heptane, a pre-exponential factor three times that of Westbrook and Dryer was 
selected. The new code uses the same kinetics (including the pre-exponential factor) as 
Westbrook and Dryer. The assumption of a single binary diffusion coefficient (along with 
neglecting pressure and thermal diffusion) results in Fick's Law for the diffusion velocity of each 
species. The binary diffusion coefficient for the fuel and oxidizer was used in the code and was 
evaluated as a function of the local temperature. Due to these assumptions, the numerical 
predictions were only in qualitative agreement with experimental results from the literature. The 
new code employs the Stefan-Maxwell equations to allow for different multicomponent diffusion 
coefficients. Employing these changes, the new code matches experimental data for the 
extinction of n-heptane droplets. 

II. Numerical Model for Droplet Combustion in a Forced Convection Environment; 

Transient Effects. 

This part of the CEFD CRT research has focused on the development of a numerical 
model to study the flame transition during fuel droplet combustion under forced convective 
environment. Based on a quasi-steady numerical model, which utilized staggered numerical grid 
and neglected heat transfer to droplet interior, a transient code has been developed to model 
droplet combustion in a forced convection, zero-gravity environment. A colocated numerical 
grid, one-step overall chemical kinetics and a single diffusion coefficient to describe the mass 
diffusion were used in the transient model. 
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Numerical results indicate that at the same ambient temperature, the Reynolds number 
determines the flame configuration that forms during droplet ignition for both suspended and 
moving droplets. In both cases, an envelope flame was formed during droplet ignition for a low 
initial Reynolds number (such as Re = 10), and a wake flame was observed at higher Reynolds 
numbers (such as Re = 50). In the suspended droplet combustion case (Figure 1), the wake flame 
configuration was observed throughout the droplet lifetime, while for the moving droplet case 
(Figure 2), as the initial wake flame gradually approached the droplet, an envelope flame might 
eventually be formed due to deceleration of the droplet. Reynolds number not only determines 
the flame configuration, but also affects the initial flame location (ignition location). 

The results for moving droplet combustion indicated that both initial Reynolds number 
and ambient temperature affect the ignition location. A high initial Reynolds number tended to 
make the droplet ignite further away from the droplet in the downstream direction. For higher 
ambient temperatures, the flame tended to initiate at locations closer to the droplet. The transient 
model predicted a marked difference between the behavior of suspended and moving droplets 
and shows the effects of both initial Reynolds number and ambient temperature on the flame 
configuration and ignition location. 

Research Outcomes 


The CEFD CRT meets weekly to allow researchers to present and discuss their new 
results. This ensures that research objectives are being met. For additional dissemination of 
findings, this CRT is in formal collaboration with both the John Glenn Research Center at Lewis 
Field in Ohio and the U.S. Department of Defense. Continued communication is a priority for 
the CEFD team. Additionally, Principal Investigator Gogos’ has maintained direct 
communication with Dr. Vedha Nayagam in San Diego, California. This communication allows 
the CEFD CRT’s modified numerical code to be validated by Dr. Nayagam’s experimental data. 
Dr. Vedha Nayagam serves as the CEFD technical monitor through his experiences at the 
National Center for Microgravity Research. The CEFD team’s numerical model is validated 
using experimental data obtained by Dr. Nayagam. 

A post doctorate research associate, Dr. Daniel Pope, and a research assistant professor, 
Dr. Hongtao Zhang, are also participants in this research project. Both researchers have set goals 
of becoming tenure track faculty in institutions of higher education. The weekly research 
meetings provide both Dr. Pope and Dr. Zhang with invaluable experience on graduate student 
advising. Pope and Zhang are strongly involved in every aspect of the CEFD research faculty 
such as, writing proposals, writing papers, presenting conference papers, and reviewing papers. 
Such mentoring opportunities are expected to continue throughout the lifetime of the research 
study. 

NASA EPSCoR funding has allowed for a variety of CEFD research activities. Those 
include the financing of travel to present papers as well as to attend conferences and workshops. 
Such papers and other substantial efforts are listed in the CEFD Team Outcomes (Appendix A). 
Student assistantships and computer services have also been made available. Additionally, the 
CEFD CRT investigations have resulted in the conference of the degree Master of Science to three 
research assistants: S.H. Soh, Y. Shi, and K. Lu. 
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Conclusion 


The Numerical Simulation of the Combustion of Fuel Droplets study is one of three 
Collaborative Research Teams (CRT) currently supported by the NSGC & EPSCoR programs. 
Each CRT strives to provide the most current information to interested members of the academic 
world. The Numerical Simulation of the Combustion of Fuel Droplets study is an evolving 
project. Periodic updates are available on a quarterly basis. 

Additional collaborations are sought with other organizations on a continual basis. All 
opportunities for collaboration are invited for consideration. Additionally, NSGC & EPSCoR 
welcome any input on program directions as well. The partnership between the NSGC & 
EPSCoR programs allows for the selection of outstanding research projects that positively 
impact aeronautical technology advancement. Those in the NSGC & EPSCoR program, the 
Collaborative Research Teams, and the industry look forward to experiencing the same high 
level of achievement in the future. A listing of CEFD CRT Team Outcomes can be found in 
Appendix A. 

Recent Activities 

1 . Weekly CRT meetings to present and discuss new results. 

2. This CRT is in formal collaboration with the John Glenn Research Center at Lewis Field 
in Ohio. 

3. This CRT is in formal collaboration with the U.S. Department of Defense. 

4. Principal Investigator Gogos’ has maintained direct communication with Dr. Vedha 
Nayagam in San Diego, California. 

5. Dr. Vedha Nayagam serves as the CEFD technical monitor through his experiences at the 
National Center for Microgravity Research. 

6. Dr. Daniel Pope and Dr. Hongtao Zhang are gaining invaluable experience on graduate 
student advising through CRT weekly meetings. 

7. This research has allowed for the financing of travel to present papers as well as to attend 
conferences and workshops. 

8. Student assistantships and computer services have also been made available through this 
CRT. 

9. The CEFD CRT investigations have resulted in the conference of the degree Master of 
Science to three research assistants: S.H. Soh, Y. Shi, and K. Lu. 
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Further Information 


If you would like further information regarding the Numerical Simulation of the 
Combustion of Fuel Droplets, or would like to collaborate with the Nebraska CRT, please 
contact: 


Dr. George Gogos, Principal Investigator 
Professor & Principal Investigator 
University of Nebraska - Lincoln 
N104 WSEC 
Lincoln, NE 68588 

TEL: 402/472-3006; FAX 402/554-1465 
E-mail: ggogos@unl.edu 


Dr. Brent Bowen, Director 
NASA Nebraska Space Grant & EPSCoR 
University of Nebraska at Omaha 
6001 Dodge Street; AH 422 
Omaha, NE 68182 

TEL: 402/554-3772; FAX 402/554-3781 
E-mail: nasa@unomaha.edu 


This overview was prepared and regularly updated by Dr. George Gogos, Dr. Brent Bowen, and 
Mrs. Jocelyn Nickerson, with contributions from various CEFD CRT members. 
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Orientation of single droplets for experiments 
Orientation of single droplets within practical combustion systems 
Numerical results for the temperature distribution surrounding a n-heptane fuel 
droplet that is moving through a combustion chamber at a velocity of 3.25 m/s 
Numerical results for the temperature distribution surrounding a n-heptane fuel 
droplet that is moving through a combustion chamber at a velocity of 5.86 m/s 
Numerical results for the temperature distribution surrounding a n-heptane fuel 
droplet that is moving through a combustion chamber at a velocity of 3.25 m/s 


36 


Figure 1 : Orientation of single droplets for experiments 
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Figure 2: Orientation of single droplets within practical combustion systems 
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Figure 3: Fuel droplet velocity 3.25 m/s 
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Figure 4: Fuel droplet velocity 5.86 m/s 
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Figure 5: Fuel droplet velocity 8.13 m/s 
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